Heat transport of clean spin-ladders coupled to phonons: 
Umklapp scattering and drag 
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We study the low-temperature heat transport in clean two-leg spin ladder compounds coupled 
to three-dimensional phonons. We argue that the very large heat conductivities observed in such 
systems can be traced back to the existence of approximate symmetries and corresponding weakly 
violated conservation laws of the effective (gapful) low-energy model, namely pseudo-momenta. 
Depending on the ratios of spin gaps and Debye energy and on the temperature, the magnetic 
contribution to the heat conductivity (K mag ) can be positive or negative, and exhibit an activated or 
anti-activated behavior. In most regimes K mag is dominated by the spin-phonon drag: the excitations 
of the two subsystems have almost the same drift velocity, and this allows for an estimate of the 
ratio Kmagj 'ttph of the magnetic and phononic contributions to the heat conductivity. 



I. INTRODUCTION 

Recent experiments both on spin-chain 1 and ladder 
compounds 2 showed a surprisingly large magnetic con- 
tribution to the heat conductivity: the heat conductiv- 
ity in the direction parallel to the ladder (attributed to 
magnons and phonons) largely exceeds the heat conduc- 
tivity in directions perpendicular to it (attributed to the 
phonons alone). The heat conductivity of clean gapless 
spin 1/2 chains coupled to phonons shows a simple ex- 
ponential behavior 3 , associated with a single character- 
istic energy scale resulting from the high-energy process 
needed to relax momentum. In contrast, gaps open up in 
the spectrum of (two-leg) spin ladders, and consequently 
the heat transport involves a complex interplay between 
different energy scales leading to a rich gamut of possible 
behaviors. In this paper, we present a theoretical frame- 
work to describe low-temperature heat transport in such 
clean gapped quasi one-dimensional systems when they 
are coupled to phonons. 

In the absence of disorder, heat transport in quasi 
one-dimensional systems is determined by momentum 
conservation (or more precisely by 'pseudo-momentum' 
conservation 3-5 ). In a clean lattice, momentum transfers 
are quantized and therefore the momentum can only de- 
cay via an Umklapp process involving a large- momentum 
high-energy state. This implies that transport in such 
systems is non-universal as it depends on both high- and 
low-energy features. Therefore controlled analytic cal- 
culations are usually not possible (the situation is, how- 
ever, simpler for systems with a finite magnetization 3 ' ). 
Nevertheless, we shall show that under certain circum- 
stances such calculations are possible. Indeed, in many 
spin-systems, the typical spin- velocity, v s , is large com- 
pared to the sound velocity of the acoustic phonons, v p . 
Therefore, the large momentum state with the lowest en- 
ergy (required for an Umklapp process) will have most of 
its momentum carried by phonons. This has two conse- 
quences: heat transport is (i) dominated by spin-phonon 



scattering and is (ii) determined by high-energy features 
of the phonon system but low-energy properties of the 
spin system. The latter observation implies that con- 
trolled calculations of transport are in principle possible 
(up to non-universal prefactors describing the electron- 
phonon coupling) using the fact that high-energy proper- 
ties of the weakly interacting phonons are often known or 
can be measured. The necessary low-energy correlators of 
the gapped spin-system can be obtained from an effective 
field theory which can be analyzed by semi-classical 6 or 
form-factor 7 methods for temperatures T below the gap. 

In contrast, in pure systems where phonons are absent 
(e.g. cold atom realizations), or when v p > v s , no such 
controlled calculation is possible as little is known about 
the non-universal high-energy properties of strongly in- 
teracting spin systems and one can make only qualita- 
tive statements, e.g. that the heat conductivity is expo- 
nentially large but finite 3 . Previous numerical studies of 
pure ladder systems at high T indicate that the heat con- 
ductivity of spin ladders is finite 8,9 but could not reach 
temperature of the order of the gap and below. On the 
analytical side, few results are available. The field the- 
oretical treatment of Ref. [10] ignored the role of Umk- 
lapp processes, thus leading to ballistic transport in a 
clean and pure system. Rozhkov and Chernyshev 11 in- 
cluded the effect of disorder and phonons in spin chains 
within a Boltzmann equation approach, but did neither 
consider spin-phonon drag nor addressed the question of 
Umklapp, which becomes essential in a clean system. 

In the following, we will first present the general field 
theoretical framework on which our calculation is based. 
We discuss the various Umklapp processes that induce 
(pseudo-) momentum decay in quanta of size G or G/2, 
where G is the reciprocal lattice vector. We shall dis- 
cuss the role of such processes in (weakly) violating the 
conservation laws of the model, allowing a hydrodynamic 
description of the system. We then use the memory ma- 
trix formalism to calculate the heat conductivity n in the 
various regimes. The memory matrix approach is gener- 
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ally not exact even in the limit of small couplings but 
it can be shown 12 to give a lower bound to k, which is 
saturated in the limit of a large separation of time scales 
between slow and fast modes, i.e. if a hydrodynamic 
description is possible, as is the case here. Finally, we 
interpret our results in terms of the spin-phonon drag 
and the decay rates of various slow modes. The various 
resulting behaviors of n{T) are summarized in Fig. 1, 
where regime I is consistent with the particular data of 
Rcf.[2] for T below the spin gap, but still larger than 
a lower scale T*, where the magnetic heat conductivity 
displays an activated behavior. 

II. LOW ENERGY EFFECTIVE THEORY 

Our starting point is the following Hamiltonian 

H = Y,H s {l) + H p + H SiP (1) 

7 

which describes an array of spin-ladders (denoted by 
's', with 7 the ladder index) coupled to acoustic three- 
dimensional phonons ('p'), via the term H sp . For sim- 
plicity, we assume H p w a u^klc^ a a k a where the 
velocities of the various branches of acoustic phonons 
(v a ) are approximated by a characteristic velocity v p , 
associated with the Debye energy via 0£> ~ v p G/2. 
A single spin ladder is described by 

Hs = J \\ E Se >i ' Se >3+ 1 + Jj - ^ Sl J ' S 2J' ( 2 ) 

where Sej is a spin- A operator acting on site j and on 
leg 1 = 1,2 of the ladder. 

As we are interested in the heat conductivity at low T, 
it is useful to consider the effective low-energy theory for 
H s (assuming a small gap, A <C J\\) described in terms 
of four massive Majorana fermions 13 

H s = I? s + 5> f dxO t (x) (3) 

In the above expression, the operators Oi are all irrel- 
evant and marginal operators allowed by the symme- 
try of the original lattice Hamiltonian (2). The three 
Majorana fields £ a , a — 1,2,3, describe the low lying 
magnon triplet with the velocity v a — v\ and (spin) gap 
A a = Ai, while the remaining Majorana field, £°, de- 
scribes a singlet excitation with gap A , and velocity 
vq: the single-particle excitations have dispersion rela- 
tion e a (k) = yjvlk 2 + A 2 . While A /Ai - 3 for weak 
J±, this ratio changes for larger J±/J\\ or when other 
microscopic interactions are present - so that we shall 
consider it as a free parameter, yet assuming Ao > Ai. 



The total heat current, obtained through the continu- 
ity equation of the energy density, is 

Je « v 2 p P p + v%P + v 2 ^ (4) 

where P p , Pq and Pi arc the momentum operators of 
phonons, singlets and triplets, respectively, for exam- 

pie Pi - -lEa=i...3, 7 /^ (^A^R + e^J- ^ 
Eq. (4) we neglect further contributions from the interac- 
tions which turn out to give only subleading corrections. 

The crucial observation on which our following analysis 
is based, is that the effective low-energy theory of our 
initial Hamiltonian conserves the total momentum, 

P T = P P + P + P 1 . (5) 

A direct consequence is that - within this low-energy 
description - the heat conductivity is infinite 14 ' 15 . How- 
ever, the continuous translational invariance of (3) is not 
a true symmetry and follows from neglecting Umklapp 
terms whose inclusion leads to a decay of P T , which now 
acquires an exponentially long life-time at low tempera- 
ture. This is to be contrasted to all other decaying modes, 
whose life-time behave as power laws of T, and allows for 
a hydrodynamic description based on the slow modes P a , 
a = p, 0, 1. Including these Umklapp operators is thus 
the correct way to obtain a low-energy effective theory 
suited to the calculation of transport properties. 

Three different Umklapp terms turn out to be impor- 
tant if we assume that v p < v , vi w : 

H u = E4 K) / dx(d x q) n co S (Gx) (6) 

n>3 J 

Hu = 9% J dx(d x q) 2 D (x)cos(Gx/2) (7) 
+ * E 9i P I dx (8 x q) 2 ^ cos(Gx) (8) 

a J 

with q(x) oc f I e ik * x (al + a_ v )-£^ the displace- 

ment field for acoustic phonons projected along the lad- 
ders (ajj is an abbreviation for the sum of contributions 
from various phonon modes). Ob = Yl a =o aa 1S tnc con " 
tinuum limit of the dimerization operator (— ) J (Si j j + 
S2j) ■ (Si j+i +S2.J+1) , where each a a (a = 0..3) is the 
Ising spin operators for the quantum Ising model that 
is naturally associated to the Majorana £ a theory 13 (in 
our conventions, Ising model '0' is in its quantum disor- 
dered phase while the three other models are ordered). 
We checked that linear couplings to phonons in (8) are 
subdominant 21 . 

On top of the Umklapp terms it is also important to 
include normal processes 

K P k = ^9a' k fdx(d x q) k QQ, (9) 

a J 

which allow momentum exchange between the spin and 
phonon systems. We do not consider normal operators 
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acting only in the spin sector or only in the phononic 
sector, as they commute separately with the spin and 
phonon momentum operators, and accordingly do not 
contribute to leading order to the heat conductivity. 



III. HYDRODYNAMIC APPROACH 

To obtain a hydrodynamic description we first iden- 
tify a basis in the space of slow modes, in our case 
given by P p ,Pq and Pi. Then we introduce the ma- 
trix of static susceptibilities of the slow modes, Xij = 
-^(Pi(O)Pj(O)) (equal time correlator), and a matrix 
of conductivities defined by Kubo formulas for the Pi's, 

= lIo /TdTeiUT ( T r P i(^ P j^)), i,j =P,0,1. The 
heat conductivity is then given by 



K — — ^ ^ Vj (Tij — K ss -\- 2/€ S p + Kp 



(10) 



Note that besides the spin and phonon heat conduc- 
tivity k ss , K pp there is also the drag term 2n ps = 

T Ei,=0,l V p V b Gbp- 

Within the memory-matrix approach 16,17 , the matrix 
of conductivities <t(uj) is expressed as 



*(w,T)=x 



1 



M(u) — iu)x 



X 



(11) 



where M(u) is the so-called memory matrix. It can be 
shown that to leading order in the coupling constants of 
the Umklapp terms, g u , and of the normal spin-phonon 
term, g N , the memory matrix is simply given by 3-5 



M, 



^((PiPjfM-iPiPjfioj), (12) 



where (...) are the retarded correlation functions evalu- 
ated with respect to the unperturbed Hamiltonian as Pj 
is already linear in the perturbations. 

To evaluate (12) to leading order, we need various 
correlation functions of the decoupled spin-phonon sys- 
tem. As discussed above and checked below, for v p < v s 
and low T, one needs high-energy correlation function 
in the phonon sector but only the low-energy asymp- 
totics of spin-correlation functions. To obtain the cor- 
rect low-energy correlators in the spin-sector it is in gen- 
eral necessary to take into account the (unitary) scatter- 
ing of the thermally excited quasi-particles using gener- 
alizations of Sachdev's semi-classical arguments 6 - see 
appendix A. For example, the correlator £/z)(x,i) = 
(Od{x, t)Oo(0, 0)) of the dimerization operator, which 
is related to the Majorana fields in a highly non linear 
and non local way, is given by: 



■i*(3x/Ci,3t/n) 



(13) 



with TV a non universal prcfactor, &(x, t) = 
x erf (x/i^/Tr) +te~ x l^ 1 \ Ko the modified Bessel func- 
tion, n = ^ e Al / T , and Ci = ViJ^KTT e Al/T - :t will, 



however, turn out that the scattering from other ther- 
mally excited quasi-particles described by Gd{x, t) is not 
important as our problem is dominated by spin-phonon 
scattering. 

A generic Umklapp memory matrix entry at zero fre- 
quency can be cast in the form g 2 9m J tdt J dxe tqx Q p Q s , 
with q = G/2 or G and Q p ( s ) the appropriate phonon 
(spin) correlator. At low T < v p q, we evaluate it in the 
saddle point approximation, by deforming the contour 
in the complex plane. The saddle point lies at one of 
the phonon propagator poles, v p t* = sign(q)x* = — 
well within the spin light-cone where there are no con- 
tributions from the spin-spin scattering. Moreover, the 
semiclassical approximation for Gd is valid at the saddle 
point. 

The memory matrix can be split into four contribu- 
tions: 



M = M N + M ph + M G/2 + M G 



(14) 



corresponding to the relaxation processes described by 
the terms (9), (6), (7) and (8) respectively. 

Generically, a memory matrix entry bears an activated 
form at low temperature, and thus can be specified by an 
activation gap and a prefactor. While the prefactor de- 
pends on non-universal, high-energy features of the spin- 
system, in the soft phonon limit v p < v s only universal 
features of the spin system determine the activation gap. 
This is why we only give the leading, exponential behav- 
ior at low T for the various memory matrix entries, (with 
the exception of one limiting case, see below). These ex- 
pressions involve different Umklapp gaps: 



E b 



v p G v p G A / 

— - E D = — + — V 



Ao 
2 



'1 



a, 



0' 



v P G 



Am/1 



(15) 



where a b 



< 1 and b = 0, 1 labeling the spin sin- 



glet and triplet excitations. These formulas have a sim- 
ple interpretation: under the constraints of energy- and 
momentum conservation they are the lowest energies for 
processes where the momentum q = G or G/2 is ab- 
sorbed by scattering a phonon from —q/2 ± 0(A/v p ) to 
q/2^f 0(A/v p ) while scattering (E a ) or creating (Ed) a 
spin excitation. The simple form of the phonon energy 
scale appearing in Eq.(15) is a result of our simplified 
treatment of the phonon sector - high energy features 
thereof would modify this form but not the fact that it 
depends of a single scale E v <~ 6d . 

The leading T-dependence of the various contribu- 
tions to the memory matrix is summarized as follows. 
The phonon Umklapp (6) gives rise to a single non van- 
ishing entry, <~ e~ Ep / T . Entries due to normal 
processes are M pv = x a + x\, Mj% = —M pb = x b , 
x b ~ e~ Ai> / T . Spin-phonon Umklapp processes with 

Vo + 2/1, 
< 1 



momentum transfer G entries are: M G p 



C b ~ 2 M G = y b with C b 



2A h 
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FIG. 1: Schematic plots of k(T) in three regimes. Depending 
on the size of E v (i.e. the Debye energy) compared to the 
values of the spin gaps one obtains qualitatively different be- 
haviors of the heat conductivity (shown for T < Ai). Note 
that in regime I and II, Kmag is displayed, while in regime III, 
the total k is displayed, (n ma g is not defined in this regime - 
see text). Below the horizontal axis are shown schematically 
the dominating momentum transfer processes within the spin 
and phonon system, and from those modes via Umklapp scat- 
tering to the lattice. T» is of order (v p /vi) 2 Ai. 

and j/6 <~ e~ Eh l T . Finally, the relevant G/2 entries are 

With all these elements at hand, using Eq. (10,11,12) 
and performing a straightforward matrix inversion yields 
the heat conductivity, which displays very different be- 
haviors depending on the relative value of the magnetic 
gaps and Dcbyc energy. In the following, we follow a more 
physically transparent approach and present in detail 
the temperature dependence of k as well as the physical 
mechanism explaining it, in different parameter regimes. 

IV. RESULTS AND DISCUSSION 

For the interpretation of our results detailed below, it 
is useful to rewrite the linear-response relation (Pj) = 
(TijAj where Aj is an external field coupling linearly to 

Pj. Using Eq. (4), we can identify it with Aj = vj^f^ 
where VTj is a fictitious temperature gradient coupling 
only to the subsystem j. Using Eq. (11) one obtains 

§ l (Pi)-xivf^ = -m- 1 h(p j ) (i6) 

where we used \ij ~ Xi^ij- 

This equation has a simple interpretation: it is a rate 
equation for the momenta and r _1 = Mj^ 1 can there- 
fore be identified with the matrix of relaxation rates. The 
matrix of conductivities can be extracted from the equi- 
librium solution of the rate equation. 

We will now discuss our results in the three different 
regimes depicted in Fig. 1, assuming always Ai < A 
and v s > v p . To investigate the relation to experiments, 
it is useful to define the magnetic contribution to k 

Kmag = K ^ph = ^ ss ^^ps ~t~ ^K pp (1^) 

where K^ h is the conductivity of a hypothetical system 
without magnetic degrees of freedom (estimated in exper- 
iments from fits to k perpendicular to the spin ladders). 



Note that there is also a negative contribution from the 
change of the phonon conductivity Sk pp = k pp — due 
to scattering from spin excitations. 



A. regime I 

We first investigate regime I of Fig. 1 where E p ,T < 
Ai and momentum relaxation is dominated by phonon 
Umklapp (E p < E\, Ed). As Ao > A l7 we can neglect the 
singlet mode and focus on the 2x2 matrix describing the 
relaxation of phonon- and triplet momentum to obtain 

_ 4^PP _ 1 4 U ^ E p /T (lg) 
K pp — Up ^ — T Xp u p T p T<Ai e ^ 10 / 

k ss = vi^- = ^ Xl vir(L P -e-^ T (19) 

n ps = vlvf-^ = f^Kpp^e^-^ T (20) 
1 Xp v p 

where x% = Xu (to leading order, \ IS diagonal) and 
1/Tp = MPp/xp ~ e~ Ev l T is the phonon Umklapp rate. 
The "pure spin heat conductivity" k ss - correspond- 
ing to heat being carried by magnons - is determined 
by the exponentially small number of spin excitations, 
Xi ~ e~ Al ' T , and the non-exponential rate of momen- 
tum transfer, l/r([^ p — M^/xi from the spin system to 
the phonon system. Most interesting is the drag term 
K psi which dominates over k ss . As we can neglect mo- 
mentum dissipation within the spin-system, even a small 
coupling of the spin and phonon systems by normal pro- 
cesses induces 'perfect drag': both subsystems have the 
same drift velocity. Therefore the ratio of heat currents 
is given by the ratio of energy densities (pf) = vfxi 22 > 
and we find 

Kps = Kpp (P E 1)/(P E P ) (21) 

in complete agreement with our memory matrix analysis. 

In Eq. (18) only the leading behavior of K pp is 
shown. To calculate Sk pp one has to consider the lead- 
ing correction arising from the subdominant mixed spin 
phonon Umklapp (8) described by the rate (t^,) -1 = 
£« M iilXv We obtain 

SKpp = -^Xpvt { 4^~^ (2Ep - El)/T ( 22 ) 

1 'sp 

This term corresponds to 'magnetic friction': the heat 
current carried by phonons is reduced due to scattering 
on the dilute gas of magnons. 

The competition between the negative Sn pp and the 
positive K ps leads to a complex cross-over behavior of 
the magnetic heat conductivity K mag , see Fig.l. Below 

T* = (2E p -E 1 )-(E p -A 1 ) w (for a x < 1), the 

negative 5n pp dominates. For Ai >7 1 >T, t , the behavior 
is fixed by the pre-exponcntial factors ; in this regime we 
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find 



^ps 



\5n p 



N_ 

on 



Ai 



T 

aT 



-9/2 



9 U ph ' 

J sp—ph 



(23) 



where the g^'s are the dimensionless microscopic cou- 
plings and TV is a numerical constant. Interestingly, the 
cross-over scale corresponds to the temperature where 
the saddle point location crosses the thermal wavelength 
of the magnons. 

Under the hypothesis that for T > T*, magnetic fric- 
tion is dominated by the enhancement of conductivity 
due to the addition of heat carriers, namely the magnons, 
we can relate the ratio of the magnetic versus phononic 
contributions to the heat conductivity, which is of experi- 
mental relevance, to purely thcrmodynamical quantities: 



^mag 

To 

K P h 



lP"Y 
^pA.pp 



(24) 



B. regime II 



Next, we consider regime II of Fig. 1, where Ai < 
E p < Aq. In this regime, the formulas (18,20,22) for 
Kpp,Kps and Supp remain unchanged, however, in con- 
trast to regime I, both n ps and Sk pp are exponentially 
large. Furthermore, the momentum transfer rate 1/t p ^\ 
from the phonon system to the spin system is now larger 
than the phonon Umklapp rate and therefore the 

two subsystems equilibrate before loosing momentum via 
phonon Umklapps, yielding: 



P 



,(B P -2A 1 )/T 



(25) 



n ss is always smaller than n pp , and generically n ss <C n ps - 
In this regime, the same remark as in regime I, regarding 
the relative size of the magnetic versus phononic heat 
conductivity, applies: if magnetic friction is subdomi- 
nant, the ratio only depends on thcrmodynamical quan- 
tities. 



C. regime III 

In regime III, E p > A , a new scattering channel (7) 
dominates where (pscudo-) momentum in quanta of G/2 
is transferred to the lattice in a complex process involv- 
ing singlet, triplet and phonon modes. The associated 
Umklapp gap Ed < E p therefore replaces E p in formulas 
(18,20,25). The leading term is the phonon contribution 



^PP ~ J, U pXp' D C 



(26) 



with = J^ij MfJ 2 jxp the rate of total momentum 

relaxation. As the phonon contribution is strongly re- 
duced by the presence of the magnetic system, the naive 



procedure of disentangling the magnon contribution n mag 
by subtracting the phonon background [Eq. (17)] is not 
useful in this regime. 



D. comparison to experiments 

The experimental data presented in Ref. [2] on cupratc 
ladders clearly display an activated behavior for K mag , 
the activation energy being close to the gap value (which 
is of the order of 400K). Measurements have been carried 
on for temperature ranging from a few tens of kelvins up 
to 300K. Due to the operative way to determine the mag- 
netic heat conductivity, namely by subtracting a 'pure 
phonon' contribution obtained by a low temperature fit, 
these data are reliable only for not too low temperatures, 
i.e. T > 50K. 

Our theory does not directly apply to these systems as 
the low-temperature properties are dominated by disor- 
der. Nevertheless, it is possible to obtain a qualitative 
understanding what happens in a situation where the 
phonon sector is disordered. Then, due to the scatter- 
ing from defects (in combination with normal phonon- 
phonon processes 18 ), the phonon heat current J^ h has 
no longer an exponentially large lifetime. It is possible 
to mimic this situation by taking the limit E p — > in 
our equations. This indicates that these compounds are 
located deep into regime I, and the corresponding acti- 
vation energy for T > is then close to the spin gap. 
Unfortunately, the temperature <~ AK turns out to be 
much lower than the minimal temperature above which 
the spin contribution to the heat conductivity can be de- 
termined reliably in the experiment. Therefore the char- 
acteristic sign change shown in the first panel of Fig. 1 is 
not observable in these samples. 



V. CONCLUSION 

Our results emphasize the richness of transport phe- 
nomena, the actual low T heat conductivity depending 
in a subtle way on the different scales present in the sys- 
tem. In particular, we find that in a clean system the 
magnetic contribution to the conductivity does not dis- 
play the trivial activated behavior with activation energy 
the spin gap. An exception is the limit E p — > which 
mimics a situation where momentum predominantly re- 
laxes via disorder in the phonon system, as might be the 
case in the cuprate systems of Ref. [2] . Wc hope there will 
be in the future more measures on spin-ladder materials 
allowing to probe the different regimes discovered by our 
study. 

The authors would like to thank S. Sachdev and A. So- 
logubcnko for stimulating discussions, and the German- 
Israeli Foundation (GIF), the DFG under SFB 608 and 
the NSF under DMR 0312495 for financial support. 
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APPENDIX A: LONG DISTANCE SPIN 
CORRELATORS 

In this appendix, we present the calculation for the low 
temperature dimerization operator correlator, Qd{x, t) — 
{Od(x, t)Oo(0, 0)), which relies on the semi-classical ap- 
proach developed in Ref. [6]. We first recall the results 
obtained by Sachdev and collaborators for a single Ising 
model, and we then apply these ideas to our system, 
which consists of four weakly coupled Ising models. 

1. Single Ising model 

The central idea of this approach is to exploit the fact 
that at low enough temperatures in a gapped system, 
r« A, typical configurations of the system correspond 
to a very dilute gas of quasiparticles, with mean quasipar- 
ticle separation of the order of e A / T , much larger than 
their thermal de Brooglie wavelength \ th = ^ V AT - As 
shown in Ref. [6] , a classical treatment of these quasipar- 
ticles is legitimate. Only the scattering rates of two par- 
ticles (diluteness allows to consider only two-body colli- 
sions) has be calculated from quantum mechanics. In one 
dimension, particles with a quadratic dispersion are per- 
fectly reflected from any potential in the limit of small 
momentum. Therefore the two-body scattering matrix 
tends to its so-called super- universal form, S = — 1, irre- 
spective of the form the two-body potential. Remarkably, 
these ingredients are sufficient to allow for a closed form 
evaluation of the Ising operator correlator, both in the 
ordered and disordered phase, with the result: 

G ord = (v(x,t)v(0,0)) = N 2 e -*(*/C.'/r) 
$(x,t) = xcrf (jj=) + te~ s '^ (Al) 

with N = (a) the local T = magnetization, ( = 

v \/jST eA ^ T ; an d T = w eA ^ T ■ I n the disordered 
phase, the correlator reads: 

G dls = (a(x, t)a(0, 0)) = A K (fVx' - vH*) e -»(*/C*/r). 

(A2) 

We are interested in the correlator of the dimeriza- 
tion operator Or> = Ila=o a<1 ■ The cru dest approx- 
imation consists of neglecting the interaction between 
the Ising models, and leads to a simple product form, 
Qd = G° dis Ila>o Gord- Of course, such an approxima- 
tion is expected to fail to capture the correct long time 
and space limit, since it neglects the scattering between 
quasiparticles on different Ising models. In reality, the 
Ising models are coupled - to lowest order, this cou- 
pling is given by the spin density-spin density coupling 
b I £r£r£l £l ~ an d the scattering is relevant in the 
sense that it qualitatively affects the form of Q D at large 
t,x/v 3> A -1 . We now proceed to take this interaction 
into account, in the limit of heavy singlet excitations, 



Aq > Ai. This is a priori the relevant physical regime, 
since in actual realizations of the spin ladder no indica- 
tion for the existence of the singlet branch, at reasonably 
low energy, has ever been reported to our knowledge. 

2. Weakly coupled Ising models 

In the semi classical approximation, field configura- 
tions contributing to the path integral representation of 
the correlation function are classical ones: quasi particles 
follow straight lines, each been given its corresponding 
Boltzmann weight. Then, the interaction between two 
particles is treated quantum mechanically, the S matrix 
bearing its super universal form (k — > limit): particles 
bounce on each other (hardcore collisions) with a ir phase 
shift (S = —1). We now repeat Sachdev et al.'s 6 line 
of reasoning including inter Ising model interaction; this 
amounts to have the straight lines representing the prop- 
agation of thermally excited states to "see" each other, 
i.e. the quasiparticles belonging to different Ising models 
to scatter, the S matrix being again the super universal 
one. 

For the Ising models n / 0, that are in the ordered 
phase, quasi particles are domain walls for the different 
Ising models that separate domains with different mag- 
netization. Each domain wall carries an index a / 0. 
Each time a domain wall a ^ crosses the line joining 
the points (0, 0) and (x, t), the operator Yl a ^ o a has its 
±1 eigenvalue flipped. This allows us to conclude that 
its semi classical correlator reads: 

( Yl ° a (x, t) J[ a a (0, 0)) = A^i 6 e -*(3x/Ci,3t/n) ( A3 ) 

It is the same as the correlator in a single ordered Ising 
model, the factor of 3 corresponding to the tripling of the 
density of excited states (the probability that a given ex- 
cited states carrying momentum p crosses the line joining 
(0, 0) to (x, t) is q' = 3q = 3\x— V\{p)t\/L, with L the sys- 
tem size and v\{p) = where ei(p) = \J A^ + v\p 2 
is the one-magnon dispersion relation). 




FIG. 2: Typical configurations contributing to Qd for differ- 
ent masses. Full lines indicate domain walls of the ordered 
models a ^ while dashed lines represent propagation of 
quasiparticle of the singlet Ising model. Left: general case. 
Right: heavy singlet limit. 

But we are rather interested in the correlator of the 
dimerization operator which includes both singlet and 
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triplet fields. The singlet Ising model being in its dis- 
ordered phase, the spin a = a 0,z has an overlap with 
the quasiparticle creation operator. Using this remark, 
the semiclassical calculation can be formulated using the 
real-space diagrams of Sachdev et al. 6 : full lines repre- 
sent world lines of domain walls in the ordered models 
; dashed lines represent the propagation of excitations 
of the disordered model ; there is a line joining (x, t) 
to (0, 0) that corresponds to the creation of a '0' parti- 
cle at point (x,t) and to its destruction at point (0,0) 
(therefore this line should be dashed on both its extrem- 
ities), and this line has to be a straight line (the state 
O D {x,t)O D {0, 0)|*) has to have a finite overlap with |<3>) 
to contribute to the trace defining the expectation value 
; therefore the pattern of world-lines with and without 
the straight line connecting (0, 0) and (x, t) has to be the 
same) - in the following we call this straight line D. Addi- 
tional rules are: each time a dashed line crosses D, there 
is a scattering event, particles bounce onto each other 
and this contributes a S = — 1 factor - note that since 
we are in one dimension and because the two scattering 
particles have the same dispersion relation, lines coming 
out of a scattering event are just continued straight. 

To proceed, one has to analyse collisions between 
dashed and full lines. A simplification occurs in the case 
of degenerate masses Ao = Ai. Then, energy-momentum 
conservation during the scattering events ensures that 
lines are also continued straight (particles just exchange 
their quantum numbers), and it can be shown that the 
problem maps onto the calculation of the Pott's spin cor- 
relator in the 4-state Pott's model in its disordered phase, 
a problem which has itself been solved in the low T regime 
using the semi-classical approach 19 . 

This fine-tuned degenerate limit is however not rele- 
vant to our situation where masses are generically differ- 
ent. If A 7^ Ai, during a scattering event between a 
domain wall and a type '0' particle, particles don't just 
exchange their momenta and lines are deviated after the 
crossing. In particular, if one considers collisions between 
domain walls and the line D, and remembering that the 
pattern of world-lines with and without the line D should 
be the same, one concludes that these events just don't 
give any contribution to the semi classical correlator - 
thus an additional rule is to forbid scattering between 
plain lines and D. Note that without D-full line cross- 
ing, the dimcrization operator is never flipped by any 
domain wall. 

Given all these rules, one has to sum up contributions 
over configurations with no crossing between full lines 
and D. In general, this is a tough problem. A simpli- 



fication occurs in the case where masses are very dif- 
ferent, and, for our purpose, when the singlet mass is 
much greater than the triplet mass. Then, it is pos- 
sible to neglect all contributions with dashed lines in 
the thermal background ; the error is of the order of 
cxp(— (Ao — Ai)/T). Enforcing that configurations con- 
tributing contain no Z?-full line crossing amounts to the 
replacement of {l-2q') N by (l-q') N in Eq.(3) of Ref.[6]. 
Performing the average over p we get: 

g D (x,t)=AK (%Vx 2 -vfr 2 ) e-**' 3 */^ 3 '/^) 

(A4) 

for x, t > 0. 




FIG. 3: Typical configuration contributing to Q' D in the heavy 
singlet mass limit. 

We will also need to evaluate correlators involv- 
ing derivatives of just one Ising spin of the kind 
{dxV a Y[ bl t a o b {x,t)0 D {Q,Q)). Of course, they arc not 
all independent. It suffices for our purpose to compute 

G' D = (a°d x l[a a (x,t)O D (0,0)), (A5) 

which we evaluate using Q' D = lim e ^o ^ Ila^o aa { x + 

e)a (x,t)O D (0,0)) - G D (x,t)j/e. Contributions to the 

first term are those with domain walls passing between 
(x+e, t) and (x, t) but not crossing the line D (sec fig. 3). 

The probability of having one single domain wall with 
momentum p passing is q e — eQ(vi(p)t — x)/L (with 
the Heaviside function), and after calculation we find 

g' D {x,t) = g D {x,t)X 1 (x,t) (A6) 

mm) = -§ r* c -«w/ r 

2 J pa n 

where po is the root of the equation dt Q^ = x/t. We 
note that |Ai(a;,i)| < fCf 1 - 
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